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Abstract 

The purpose of this article is to present, in a simple way, an analytic 
approach to special numbers and polynomials. The approach is based on 
the derivative polynomials. The paper is, to some extent, a review article, 
although it contains some new elements. In particular, it seems that some 
integral representations for Bernoulli numbers and Bernoulli polynomials 
can be seen as new. 
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1 Introduction 

Let u = u(z) be a holomorphic function defined in a domain £)„ C C which 
fulfills the Riccati differential equation with constant coefficients 

u' = r{u — a){u — b) (1) 

where r,a,b are real or complex numbers r ^ 0, a ^ b. Let v = v{z) be a 
holomorphic function defined in a domain Dy C C which is related with u(z) 
and fulfills the following differential equation 



rylu-^], (2) 



where a, b, r, u{z) are as in JT]). 

Examples of such pairs of functions and equations are: 
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1. u{z 

2. u{z 

3. u{z 

4. u{z 

5. u{z 



= tanz, = + 1, u(z)=secz, v' — vu, 

= tanhz, = — + 1, v{z) = 1/coshz, w' = - 

= cotz, u'{z) = —u^ — 1, v{z) = cscz, v' — ~vu, 
— cothz, x' {z) = — 1, v{z) ~ 1/sinhz, v' — vu, 
= l/(l + e^), u'{z) = u^-u, = 6^/2/(1 + 6^), 



v' = v{u - 1/2), 

6. = 1/(1 + e-^), u'(z) = + w(z) = 6-^/2/(1 + e-^), 
i;' ^ -~v{u - 1/2), 

7. more generally the logistic function: u{z) = q/{l +pe~*^). 



qe 



-sz/2i^l + pe"^^), ?;'(2) = -v{q/2 - u) 



(with p > 1, q > 0, s > 0). 
We will consider also the following generalization of equation ([2]) 

a + b 



(3) 



where d is a real or complex number. 

Such system of differential equations has been investigated by Hoffman [12] 
(instead of equations ([2][3]) he regarded v' = vu but in each particular case, 
which he regarded, a + b — 0) and by Franssens [10] (who investigated the 
equation v' = —vu) 

Let {fli, a2, . . . , a„} be a permutation of the set {1, 2, . . . , n}. Then {oj, flj+i} is 



The Eulerian number 



is defined 



an ascent of the permutation if aj < aj^i. 

as the number of permutations of the set {l,2,...,n} having k permutation 
ascents (see [11], p. 267). For example for n = 3 the permutation {1,2,3} has 
two ascents, namely {1,2} and {2,3}, and {3,2,1} has no ascents. Each of 

/3\ 



the other four permutations of the set has exactly one ascent. Thus 



1, 



= 4, and 



= 1. It is well known that Eulerian numbers satisfy the 



following relations: 



n 
k 

n + 1 
k 

' n 
k 



n 

n — k — 1 
n 
k 



(fc + 1) 



k 



3=0 



+ (n 
+ 1 

i 



n 
k - 1 



{k-j + ir. 



(4) 
(5) 
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The Eulerian polynomial En{x), n = 0, 1, 2, . . . is defined (see Comtet p]) by 
the formula 

There is a slightly different definition of the Eulerian polynomial A„ (x) (see for 
example Foata '8 ) i.e., 

^n(^)=E(fc)^'' ^0(2;) -1. (7) 
k=Q ^ ' 

Thus En{x) = xA„(2;) for n>l,Ai{x)^ Ao{x) ^ Eo{x) = 1. 

The MacMahon numbers {Mn k} are defined by the recurrence formula (see 

M) 

MnM = (2fc - l)M„_i,fe + (2n - 2A: + l)M„_i,fc_i, (8) 
where 1 < fc < 71, Af(n, 1) = 1, n = 1, 2, . . .. 

The MacMahon polynomial M„(a;), ri = 0, 1, 2, . . . is defined as follows 

n+l 

M„(x) = EM„+i,fea;'=-i. (9) 
fc=i 

2 Derivative polynomials 

The following theorem has been discussed during the Conference ICNAAM 2006 
(September 2006) in Greece and it appears in the paper 15 . For convenience 
of the reader we give it with an inductive proof. Independently the theorem 
has been considered and proved, with a proof based on generating functions, by 
Franssens [10] (see also [16]). 

Theorem 1. // a function u[z) satisfies equation (Qp, then the nth derivative 
of u{z) can be expressed by the following formula: 

u(")(^) = r"E/^\(M-a)'=+i(u-6)"-'^ (10) 
fe=0 ^ ' 

where n — 2,3, . . .. 
Proof. By ([T]) we get 

u"(t) = r[{u -a) + {u- b)]u'{z) = r^[{u - a)(u - b)^ + {u - a)^{u - b)], 

which establishes PH)) for n — 2. Let us assume that for an integer n > 2 
formula (fTI?|) holds. Using recurrence formula Q in the last step of the following 
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calculation we get 



, n — 1 

a \ ^ /n 



k=0 



71 — k 



71-1 



x(ii-a)'=+i(?i 



n—1 



(u-a)"+^(u-6) 



n + 1 



{u-af+'{u-b) 



fc=0 

which ends the proof. 



□ 



The following two theorems (2 and 3) are connected with solutions of equa- 
tions ([2]) and ([3]) respectively and are due to Franssens [TO]- Theorem 2 is a 
particular case of Theorem 3. We write them down here in a slightly different 
form than in |10| . Franssens proved the theorems by using generating functions 
but they can be proved also by induction, similarly as Theorem 1. 

Theorem 2. // functions u — u{z) and v — v{z) are any solutions of the 
equations (QJ), (0) respectively, then the nth derivative of v{z) is equal: 



n+l 



(11) 



fc=i 



Denote by Qn{u; a, b), n — 0,1,2, . . . the polynomial (of order n) standing 
on the right hand side of equation ([TT|) i.e., 



n+l 



,{u;a, b) = J2 Mn+iMu - ar+'~\u - bf-\ 



fc=i 



Theorem 3. // functions u = u[z) and v = v{z) are any solutions of the 
equations (QP, @) respectively, then the nth derivative of v{z) is equal: 



k=0 



(12) 



The polynomials Qn{u;a,b) are related to MacMahon polynomials © by 
the formula 

(u — oj" 
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Similarly we denote by P„+i(w; a, b), n — 1,2, . . . the polynomial (of order n + 1) 
standing of the right hand side of equation (|T0| . Thus 



Pn+i{u;a,b) = (^'^yu - a)'=+i(u ~ 6)"-^ n = 1,2,. . . = - a. 

Obviously the polynomial Pn+i (u; a, b) can be rearranged into the Eulerian poly- 
nomial En{x), n = 1,2,... using the following formula: 

Pn+i{u;a,b) 

Polynomials {P„(w; a, b)} and a, 6)}t are called the derivative polynomi- 

als. They have been introduced by Hoffman [12] who used them to calculate 
some integrals with parameters and for summing some series, without giving any 
explicit formula for the coefficients. The polynomials were recently intensively 
studied (see for example [T]. [2]. [6]. [9]. [T6 ] . [T7 ] ). 

3 Generating functions for the Eulerian polyno- 
mials 

It is easy to find the closed form of the following exponential generating function 
(see [12], [10]): 

^2 

F{u,t)=u + rP2iu;a,b)t + r^P3{u;a,b)— + ---. (14) 

For convenience of the reader we give the calculation for (fT4)) . Let u = u{z) be 
a solution of the equation ([Ij. By the Taylor formula for the function u — u{z) 
we have 

^2 

F{u{z),t) = u{z) + rP2{u{z);a, b)t + P3{u{z); a, b)— + ■■■ 
= u{z) + u'{z)t + u"{z)— + • • ■ = u{z + t), 

and 

F{u,t) =u{z{u)+t). (15) 

For example ifa^O, b = 1, u = l/(l + exp(2;)), exp(z) = (1 — u)/u (see point 
[5j on the list of functions on page 2) we get 

1 1 u 

^> = l + e-(u)+t = 1 + l^e* " u+il~u)e'- 

The generating function (|16p can be used for calculation of the exponential 
generating function for the Eulerian polynomials ^ 

2 3 

Eo{x) + E,{x)y + E2{x)^ + Esix)^ + ■■■■ (17) 
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In order to do it let us observe, that wc obtain the generating function (IT71) by 
substituting in the expression 

Fiu,t)-u 
u-1 ^ ' 

where F{u, t) is given by the formula (I16|) . u/ {u — 1) — x and {u — \)t = y (that 
is 1/(m — 1) = a; — 1, xl{x — 1) — li, i — (x — l)y). We calculate 



F{u,t) - u F{u,t) - 1 _ u+{i-u)e* 



u — 1 u ~ I u — 1 u + {1 — u)e* 



x — 1 x — 1 

Therefore the generating function for the Eulerian polynomials is 

Eo{x) + E,{x)y + E,ixf^ + E^ixf^ + ■■■= (18) 

Formula (|18l) gives immediately the generating function for the Eulerian poly- 
nomials {An{x)} defined by ([7]). We have 

2 3 

Ao{x) + Ai{x)y + A2{x)^ + A^ix)^ + ■■■ 

1 — x \ 1 ^ x — 1 ,^ , 

^1 - .'Ee(i-^)!' Jx ~a;-e(^-i)!' 

Foata [8] notices that formula (|19p was known to Euler. 

4 Some others classical formulae concerning Eu- 
lerian polynomials 

The approach to Eulerian numbers and polynomials presented here is useful 
in obtaining other known results. For example the following classical formula 
concerning Eulerian polynomials 

En{x) = r )E,{x){x - 1)"-!-'-^ + E,{x){x - l)"-i (20) 

k=l ^ ^ 

n = 1, 2, ... or equivalently expressed in terms of the polynomials{A„(a;)}: 

n-l , ^ 

An{x)^Y.[l)Mx){x-ir-'-' (21) 

is an easy consequence of the following lemma. 
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Lemma 1. // a function u = f{z) fulfills the equation u' = u{u — 1) then for 
any n = 1,2,... 



/(")(z) = (/(z)-i)5:r)/«(z). 

I — n V / 



(22) 



k=0 

Proof. The proof is by induction with respect to n. For n = 1 formula 
is obviously true and let us suppose that it holds for a positive integer n. We 
have 

n-l y ^ n-1 , ^ 

(z) = /'WE (fcj/^'H^) + (/w - i)E ( J/^'^'^w 

fe=0 ^ ' k=0 ^ ^ 

= (/(z)-i)/(z)Er /('=n-) + (./w-i)Eu"i (23) 

/£ = ^ ^ fe=l ^ 

By rearranging (f22|) to the form 

fc=0 ^ ^ fe=0 ^ ^ 

and using it to the first sum of (|23)) we get 



=(/w-i)E("r)^^'^(^)' 

and then formula (1^^ is proved. □ 
By using Theorem 1 we see that formula is equivalent to 

E {"^y-^u - 1)"- = (. - 1) ||e (I) E ( ■)"^"^(" - 1)"^ + - 

By substituting here — 1) = x, w = x/{x — 1), u — 1 = — 1) we get 

n — 1 / V ^-ii ^ / 71 — 1 / \ k—1 



hence we obtain the formula 
and the formulae ([201) and (|21l) are proved. 
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5 Generating functions for the MacMahon poly- 
nomials 



It is useful to get the generating function for the polynomials Q„(u) as 

G{u,t) = Qo{u;a,b) + -Qi{u:a,b)t+ ■^Q2{u;a,b)— + -^Qsiu; a,b)— H 

Let the functions u = u{z), v = v{z) fulfill respectively equations ^ and ([2]). 
Then using (fTTj) we have 

v{z)G{u{z),t) = v{z) + v{z)-Qi{u{z)]a,b)t + v{z)—Q2{v{z)]a,h)—A 

= v{z) + v'{z)t + v"{z)j^_ +---^g{z + t). 
For example, for data from the point 5 on page 2 we get 

G{u{z\t) = -—^^, 
1 + 1 + e'^e^ 

hence 

1 + e^e' 

Therefore in this case 



(1+ i_pi)e*/2 ^ et/2 
1 + i^e* ^ u+ (1 -M)e* 



, rJ . , ...r (24) 



The generating function can be used for calculation of the exponential 
generating function for the MacMahon polynomials ([9]) 

M^{x) + Uh{x)v + ^M2ix)^ + ^M:,ix)^ + ■■■. (25) 

In order to do it let us observe, that we obtain the generating function (|25|) by 
substituting into G{u,t) given by (|24|) . u/{u — 1) = x and {u — l)t = y (that is 
l/(u — 1) — X — 1, x/{x — 1) = u, t = {x — l)y). We calculate 



Thus 

Mo(x) + -M^ix)y + ^M2ix)l^ + ^M3(x)|j + • • • = \^J^^,_,^^/, 

and 

y2 3 (l-x)e'-^~^^y 

Mo{x) + Ah{x)y + M,{xf- + M^ixf- + ■■■= \_J^,,_^^y ■ 
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6 Integral representations 



In paper [17| we have proved that for n — 1,2,3,... 

/ Pn{u;a,b)du = -{b-a^+^Bn, (26) 

J a 

where -B„ is the nth BernouUi number and since Pn{u] a, b) is a polynomial i.e., 
an entire function, the integral can be understood as integral over any curve 
(piecewise smooth), joining points a and b. Formula (j26l) is important because 
it gives immediately the following Grosset-Veselov formula (see Grosset-Veselov 

(_-\\m—l r + oo / jm — 1 i \^ 

B2m = ^ / -7 T dx, (27) 

which connects one-soliton solution of the KdV equation with Bernoulli num- 
bers. Fairlie and Veselov proved, by using the conservation laws, that KdV 
equation is directly related to the Faulhaber polynomials and the Bernoulli poly- 
nomials (see [13]). Grosset and Veselov [4] demonstrated the formula (|27l) in two 
ways, using the cited results and then adapting an idea due to Logan described 
in the book [11]. Boyadzhiev [3] gave an alternative proof of ([ST]) , based on the 
Fourier transform. He noted that this proof was independently suggested by 
Professor A. Staruszkiewicz (see also 'Note added in Proofs' at the end of [4]). 

In order to prove (|27|) let us observe that one of the solutions of the equation 
([!])■ for a — —1, b = l,r = —1, is u = tanhz. Since the image of the real line, 
under this function, is the interval (—1,1) we have by ((26|) 

r ^^^"^"-^'"''^ dz = (-1)"-! Pn{u; -1, l)du = (-l)"2"+ii?„. (28) 
J-oo cosh z J-i 

Taking in p8|) n = 2m and using, (m — l)-times, the formula of integration by 
parts for the leftmost integral, we get the Grosset-Veselov formula (|27l) . 
There arises a natural question about similar calculations for other polynomials 
e.g., for Qn{u, a,b). As we have proved in Sec. 4 (formula ([24])) the generating 
function for the polynomials in the case of a = 0, 6 = 1, r = 1 is 

G{u,t)^Qo{u;0,l)) + lQi{u;0,l)t+^Q2{u;0,l)f 



2! w+(l-it)e* 



and therefore 

r-l 



However since the generating function for the Bernoulli polynomials 
Bo{w),Bi{w),. . . is 

Bo{w) + Bi{w)t + B2{w)- + B3{w)- + --- = (30) 
then from (|29)) we get the following theorem. 
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Theorem 4. For n = 0, 1, 2, 



^ g„(w;0,l)du = 2"B„ . (31) 

Since the polynomial Qn{u\a,h) is homogenous then by a suitable Unear 
change of the variable in the integral (I5T|) we get immediately 

^ g„(u;a,6)d7. = 2"S„ (&-a)«+i. (32) 

Let us recall that _B„ = -Bn(O). Then in view of (I^Bl) and p2p the next natural 
question arises, which concerns the existence of a family of polynomials 'con- 
necting' polynomials Pn{u] a, b) and Qn{u] a, b) in the sense that the correspond- 
ing integrals would give the values of the Bernoulli polynomial at intermediate 
points between and ^. 

Denote by S'„(u; a, b, d), n ~ 0,1,2, .. . the polynomial (of order n) standing on 
the right hand side of the equation i.e., 

Sn{u; a, ^' = Q {'^d)''Qn-k{u; a, b). (33) 

We will prove that {S'„(u; a, b, d)} form the requested family of polynomials. A 
closed form formula for the following exponential generating function: 

H{u,t) = Sq{u; a, b,d) + -Si{u;a, b,d)t + -^S2{u; a, b,d)-^^ H , (34) 

can be found similarly as in the previous cases. We assume that functions 
u = u{z) and v ~ v{z) are solutions of the equations (HJ and respectively. 
Using the Taylor formula for the function v = v{z) we have 

r 

v{z)H{u{z),t) — v{z)So{u{z);a,b,d) + v{z)-Si{u{z);a,b,d)t 
+v{z)—S2iu{z);a,b, d)^ --- 
= v{z) + v'{z)t + v"{z)- + ---=v{z + t). (35) 

For example taking here a = 0,&=l,r = l and u{z) ~ the second 



equation ([3]) has the form 
with a solution 

,(l/2+d)z 



1 + e'- 



v{z) = — . (36) 
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Using (1551) and we get 

v{z + t) (l + e^)e(i/2+<i)« 



H{u{z),t) 



v{z) 1 + e^e* 



and putting here u(z) — , e — we arrive at 

1 + u 



Then dSZl) yields 



(y, t) ^ ^ ^ " = _ . (37) 

^ ' ^ 1 + i^e* M+ 1-M)e* ^ 



1 /•! g(l/2+d)t ^g(l/2+ti)t 

and therefore using (pO| we arrive at the formula [n = 0,1,2,.. .) 



(38) 



In order to generahze ([55)1 to the polynomial Sniu;a,b,d) we use formula ([55]). 
Therefore by ([5^ we have 

J Sniu;a,b,d) ^J2(^{2d)'' Qn-kiu; a,b)du 

k-- 



b ~ a 
d 



2 6 



At the very end of the above calculation we have used the addition formula for 
the Bernoulli polynomials (see Temme [18], p. 4) 



k=0 

Thus we have proved the following 
Theorem 5. For n = 1,2, . . . 

rb 



^ 5„(u;a,6,d)dM = 2"(6-a)"+iB„ Q 



b — a 
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Comparing the generating functions (with parameters a ~ 0, b = 1, r = 
1, d = -1/2): F{u, t) (given by formula ^) with H{u, t) (formula dST])) of the 
polynomials {P„(w; 0, 1)} and {Sn{u] 0, 1, —1/2)} respectively we get also 

F„+i(.;0,l) ^^ 

u 2" ' ' ' ' ^ 

In particular, it follows that the coefficients of the polynomial 
— Sn{u\ 0, 1, —1/2) are all integer. 
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